Concurrent impact of de novo mutations on cranial and cortical development in nonsyndromic craniosynostosis.

OBJECTIVE: Nonsyndromic craniosynostosis (nsCS), characterized by premature cranial suture fusion, is considered a primary skull disorder in which impact on neurodevelopment, if present, results from the mechanical hindrance of brain growth. Despite surgical repair of the cranial defect, neurocognitive deficits persist in nearly half of affected children. Therefore, the authors performed a functional genomics analysis of nsCS to determine when, where, and in what cell types nsCS-associated genes converge during development. METHODS: The authors integrated whole-exome sequencing data from 291 nsCS proband-parent trios with 29,803 single-cell transcriptomes of the prenatal and postnatal neurocranial complex to inform when, where, and in what cell types nsCS-mutated genes might exert their pathophysiological effects. RESULTS: The authors found that nsCS-mutated genes converged in cranial osteoprogenitors and pial fibroblasts and their transcriptional networks that regulate both skull ossification and cerebral neurogenesis. Nonsyndromic CS-mutated genes also converged in inhibitory neurons and gene coexpression modules that overlapped with autism and other developmental disorders. Ligand-receptor cell-cell communication analysis uncovered crosstalk between suture osteoblasts and neurons via the nsCS-associated BMP, FGF, and noncanonical WNT signaling pathways. CONCLUSIONS: These data implicate a concurrent impact of nsCS-associated de novo mutations on cranial morphogenesis and cortical development via cell- and non-cell-autonomous mechanisms in a developmental nexus of fetal osteoblasts, pial fibroblasts, and neurons. These results suggest that neurodevelopmental outcomes in nsCS patients may be driven more by mutational status than surgical technique.

Pathogenic variants in the paired-related homeobox 1 gene (PRRX1) cause craniosynostosis with incomplete penetrance.

PURPOSE: Studies have previously implicated PRRX1 in craniofacial development, including demonstration of murine Prrx1 expression in the preosteogenic cells of the cranial sutures. We investigated the role of heterozygous missense and loss-of-function (LoF) variants in PRRX1 associated with craniosynostosis. METHODS: Trio-based genome, exome, or targeted sequencing were used to screen PRRX1 in patients with craniosynostosis; immunofluorescence analyses were used to assess nuclear localization of wild-type and mutant proteins. RESULTS: Genome sequencing identified 2 of 9 sporadically affected individuals with syndromic/multisuture craniosynostosis, who were heterozygous for rare/undescribed variants in PRRX1. Exome or targeted sequencing of PRRX1 revealed a further 9 of 1449 patients with craniosynostosis harboring deletions or rare heterozygous variants within the homeodomain. By collaboration, 7 additional individuals (4 families) were identified with putatively pathogenic PRRX1 variants. Immunofluorescence analyses showed that missense variants within the PRRX1 homeodomain cause abnormal nuclear localization. Of patients with variants considered likely pathogenic, bicoronal or other multisuture synostosis was present in 11 of 17 cases (65%). Pathogenic variants were inherited from unaffected relatives in many instances, yielding a 12.5% penetrance estimate for craniosynostosis. CONCLUSION: This work supports a key role for PRRX1 in cranial suture development and shows that haploinsufficiency of PRRX1 is a relatively frequent cause of craniosynostosis.

De novo variants implicate chromatin modification, transcriptional regulation, and retinoic acid signaling in syndromic craniosynostosis.

Craniosynostosis (CS) is the most common congenital cranial anomaly. Several Mendelian forms of syndromic CS are well described, but a genetic etiology remains elusive in a substantial fraction of probands. Analysis of exome sequence data from 526 proband-parent trios with syndromic CS identified a marked excess (observed 98, expected 33, p = 4.83 × 10-20) of damaging de novo variants (DNVs) in genes highly intolerant to loss-of-function variation (probability of LoF intolerance > 0.9). 30 probands harbored damaging DNVs in 21 genes that were not previously implicated in CS but are involved in chromatin modification and remodeling (4.7-fold enrichment, p = 1.1 × 10-11). 17 genes had multiple damaging DNVs, and 13 genes (CDK13, NFIX, ADNP, KMT5B, SON, ARID1B, CASK, CHD7, MED13L, PSMD12, POLR2A, CHD3, and SETBP1) surpassed thresholds for genome-wide significance. A recurrent gain-of-function DNV in the retinoic acid receptor alpha (RARA; c.865G>A [p.Gly289Arg]) was identified in two probands with similar CS phenotypes. CS risk genes overlap with those identified for autism and other neurodevelopmental disorders, are highly expressed in cranial neural crest cells, and converge in networks that regulate chromatin modification, gene transcription, and osteoblast differentiation. Our results identify several CS loci and have major implications for genetic testing and counseling.

Multiomic analyses implicate a neurodevelopmental program in the pathogenesis of cerebral arachnoid cysts.

Cerebral arachnoid cysts (ACs) are one of the most common and poorly understood types of developmental brain lesion. To begin to elucidate AC pathogenesis, we performed an integrated analysis of 617 patient-parent (trio) exomes, 152,898 human brain and mouse meningeal single-cell RNA sequencing transcriptomes and natural language processing data of patient medical records. We found that damaging de novo variants (DNVs) were highly enriched in patients with ACs compared with healthy individuals (P = 1.57 × 10-33). Seven genes harbored an exome-wide significant DNV burden. AC-associated genes were enriched for chromatin modifiers and converged in midgestational transcription networks essential for neural and meningeal development. Unsupervised clustering of patient phenotypes identified four AC subtypes and clinical severity correlated with the presence of a damaging DNV. These data provide insights into the coordinated regulation of brain and meningeal development and implicate epigenomic dysregulation due to DNVs in AC pathogenesis. Our results provide a preliminary indication that, in the appropriate clinical context, ACs may be considered radiographic harbingers of neurodevelopmental pathology warranting genetic testing and neurobehavioral follow-up. These data highlight the utility of a systems-level, multiomics approach to elucidate sporadic structural brain disease.

Molecular Scalpels: The Future of Pediatric Craniofacial Surgery?

SUMMARY: CRISPR-Cas genome editing tools are among the most substantial advances in the life sciences in modern history. Single-dose gene therapies to correct pathogenic mutations have moved quickly from bench to bedside, with several therapeutics designed by CRISPR pioneers entering various stages of clinical investigation. Applications of these genetic technologies are poised to reshape the practice of both medicine and surgery. Many of the most morbid conditions treated by craniofacial surgeons are syndromic craniosynostoses caused by mutations in fibroblast growth factor receptor genes, including Apert, Pfeiffer, Crouzon, and Muenke syndromes. The fact that pathogenic mutations in these genes are recurrent in the majority of affected families presents a unique opportunity to develop "off-the-shelf" gene editing therapies to correct these mutations in affected children. The therapeutic potential of these interventions could reshape pediatric craniofacial surgery, potentially first eliminating the need for midface advancement procedures in affected children.

Quantifying the Impact of Genetics on Neurocognition in Nonsyndromic Sagittal Craniosynostosis.

BACKGROUND: Previous work has identified an association between de novo and transmitted loss-of-function mutations in genes under high evolutionary constraint with neurodevelopmental delays in nonsyndromic craniosynostosis (NSC). The authors sought to quantify the neurocognitive effect of these genetic lesions. METHODS: In a prospective, double-blinded cohort study, demographic surveys and neurocognitive tests were administered to patients recruited from a national sample of children with sagittal NSC. Scores for academic achievement, Full-Scale Intelligence Quotient (FSIQ), and visuomotor skills were directly compared between patients with and without damaging mutations in genes with a high probability of loss of function intolerance using two-tailed t tests. Analysis of covariance was also used to compare test scores while controlling for surgery type, age at surgery, and sociodemographic risk. RESULTS: Fifty-six patients completed neurocognitive testing, 18 of whom had a mutation in a highly constrained gene. There was no significant difference between groups in any sociodemographic factors. After controlling for patient factors, patients with high-risk mutations had poorer performance compared with patients without high-risk mutations in every testing category, with significant differences in FSIQ (102.9 ± 11.4 versus 110.1 ± 11.3; P = 0.033) and visuomotor integration (100.0 ± 11.9 versus 105.2 ± 9.5; P = 0.003). There were no significant differences in neurocognitive outcome when stratifying groups based on type of surgery or age at time of surgery. CONCLUSIONS: Even after controlling for exogenous factors, the presence of mutations in high-risk genes led to poorer neurocognitive outcomes. High-risk genotypes may predispose individuals with NSC to deficits, particularly in FSIQ and visuomotor integration. CLINICAL QUESTION/LEVEL OF EVIDENCE: Risk, II.

Damaging variants in FOXI3 cause microtia and craniofacial microsomia.

PURPOSE: Craniofacial microsomia (CFM) represents a spectrum of craniofacial malformations, ranging from isolated microtia with or without aural atresia to underdevelopment of the mandible, maxilla, orbit, facial soft tissue, and/or facial nerve. The genetic causes of CFM remain largely unknown. METHODS: We performed genome sequencing and linkage analysis in patients and families with microtia and CFM of unknown genetic etiology. The functional consequences of damaging missense variants were evaluated through expression of wild-type and mutant proteins in vitro. RESULTS: We studied a 5-generation kindred with microtia, identifying a missense variant in FOXI3 (p.Arg236Trp) as the cause of disease (logarithm of the odds = 3.33). We subsequently identified 6 individuals from 3 additional kindreds with microtia-CFM spectrum phenotypes harboring damaging variants in FOXI3, a regulator of ectodermal and neural crest development. Missense variants in the nuclear localization sequence were identified in cases with isolated microtia with aural atresia and found to affect subcellular localization of FOXI3. Loss of function variants were found in patients with microtia and mandibular hypoplasia (CFM), suggesting dosage sensitivity of FOXI3. CONCLUSION: Damaging variants in FOXI3 are the second most frequent genetic cause of CFM, causing 1% of all cases, including 13% of familial cases in our cohort.

Molecular genetics of human developmental neurocranial anomalies: towards "precision surgery".

Recent trio-based whole-exome sequencing studies of congenital hydrocephalus and nonsyndromic craniosynostosis have identified multiple novel disease genes that have illuminated the pathogenesis of these disorders and shed new insight into the genetic regulation of human brain and skull development. Continued study of these and other historically understudied developmental anomalies has the potential to replace the current antiquated, anatomically based disease classification systems with a molecular nomenclature that may increase precision for genetic counseling, prognostication, and surgical treatment stratification-including when not to operate. Data will also inform future clinical trials, catalyze the development of targeted therapies, and generate infrastructure and publicly available data sets relevant for other related nonsurgical neurodevelopmental and neuropsychiatric diseases.

De novo mutations in the BMP signaling pathway in lambdoid craniosynostosis.

Lambdoid craniosynostosis (CS) is a congenital anomaly resulting from premature fusion of the cranial suture between the parietal and occipital bones. Predominantly sporadic, it is the rarest form of CS and its genetic etiology is largely unexplored. Exome sequencing of 25 kindreds, including 18 parent-offspring trios with sporadic lambdoid CS, revealed a marked excess of damaging (predominantly missense) de novo mutations that account for ~ 40% of sporadic cases. These mutations clustered in the BMP signaling cascade (P = 1.6 × 10-7), including mutations in genes encoding BMP receptors (ACVRL1 and ACVR2A), transcription factors (SOX11, FOXO1) and a transcriptional co-repressor (IFRD1), none of which have been implicated in other forms of CS. These missense mutations are at residues critical for substrate or target sequence recognition and many are inferred to cause genetic gain-of-function. Additionally, mutations in transcription factor NFIX were implicated in syndromic craniosynostosis affecting diverse sutures. Single cell RNA sequencing analysis of the mouse lambdoid suture identified enrichment of mutations in osteoblast precursors (P = 1.6 × 10-6), implicating perturbations in the balance between proliferation and differentiation of osteoprogenitor cells in lambdoid CS. The results contribute to the growing knowledge of the genetics of CS, have implications for genetic counseling, and further elucidate the molecular etiology of premature suture fusion.

Impaired neurogenesis alters brain biomechanics in a neuroprogenitor-based genetic subtype of congenital hydrocephalus.

Hydrocephalus, characterized by cerebral ventricular dilatation, is routinely attributed to primary defects in cerebrospinal fluid (CSF) homeostasis. This fosters CSF shunting as the leading reason for brain surgery in children despite considerable disease heterogeneity. In this study, by integrating human brain transcriptomics with whole-exome sequencing of 483 patients with congenital hydrocephalus (CH), we found convergence of CH risk genes in embryonic neuroepithelial stem cells. Of all CH risk genes, TRIM71/lin-41 harbors the most de novo mutations and is most specifically expressed in neuroepithelial cells. Mice harboring neuroepithelial cell-specific Trim71 deletion or CH-specific Trim71 mutation exhibit prenatal hydrocephalus. CH mutations disrupt TRIM71 binding to its RNA targets, causing premature neuroepithelial cell differentiation and reduced neurogenesis. Cortical hypoplasia leads to a hypercompliant cortex and secondary ventricular enlargement without primary defects in CSF circulation. These data highlight the importance of precisely regulated neuroepithelial cell fate for normal brain-CSF biomechanics and support a clinically relevant neuroprogenitor-based paradigm of CH.

Virtual Neurocognitive Testing in Craniofacial Surgery.

SUMMARY: Digital communication platforms have had an impact on plastic surgery research. Although neurocognitive testing is an important tool for understanding the long-term developmental outcomes of congenital craniofacial conditions, testing has been geographically limited because of its interactive nature. Virtual neurocognitive testing is a facsimile version of in-person testing that allows researchers to overcome this limitation with high fidelity. Ten school-age subjects completed a neurocognitive assessment battery measuring intelligence quotient, academic achievement, and visuomotor skills at an average age of 7.8 ± 1.4 years. The subjects then retook the same battery through an identical virtual testing method approximately 1 year after the initial test. There were no significant differences in verbal intelligence quotient (p = 0.878), performance intelligence quotient (p = 0.813), and full-scale intelligence quotient (p = 0.982) scores obtained through in-person and virtual tests. There were no significant differences between academic achievement letter and word recognition (p = 0.999), math computation (p = 0.619), reading comprehension (p = 0.963), and spelling (p = 0.344) scores. There were no significant differences in visuomotor integration (p = 0.158), visual perception (p = 0.798), and motor coordination (p = 0.796) scores obtained through the two methods. The median at-home time to set up study equipment for the virtual testing was 5.0 minutes (interquartile range, 2.3 minutes). Seventy percent of participants reported that they preferred participating in research from home. All participants stated that the virtual research study was an overall positive experience. Virtual neurocognitive testing is an effective method of performing high-fidelity neurocognitive assessments while overcoming geographic barriers.

Genetic Influence on Neurodevelopment in Nonsyndromic Craniosynostosis.

BACKGROUND: Nonsyndromic craniosynostosis is one of the most common anomalies treated by craniofacial surgeons. Despite optimal surgical management, nearly half of affected children have subtle neurocognitive deficits. Whereas timing and type of surgical intervention have been studied, the possibility of genetic influence on neurodevelopment in nonsyndromic craniosynostosis patients remains unexplored. METHODS: The authors performed whole-exome sequencing for 404 case-parent trios with sporadic nonsyndromic craniosynostosis. Statistical analyses were performed to assess the burden of de novo mutations in cases compared to both expectation and 1789 healthy control trios. Individuals with and without each mutation class were analyzed, and the presence or absence of various types of neurodevelopmental delay were recorded alongside demographic information. RESULTS: The authors identified a highly significant burden of damaging de novo mutations in mutation-intolerant [probability of loss of function intolerance (pLI) >0.9] genes in nonsyndromic craniosynostosis probands (p = 5.9 × 10-6). Children with these mutations had a two-fold higher incidence of neurodevelopmental delay (p = 0.001) and a more than 20-fold greater incidence of intellectual disability (p = 7.2 × 10-7), and were 3.6-fold more likely to have delays that persisted past 5 years of age (p = 4.4 × 10-4) in comparison with children with nonsyndromic craniosynostosis without these mutations. Transmitted loss of function mutations in high-pLI genes also conferred a 1.9-fold greater risk of neurodevelopmental delay (p = 4.5 ×10-4). CONCLUSIONS: These findings implicate genetic lesions concurrently impacting neurodevelopment and cranial morphogenesis in the pathoetiology of nonsyndromic craniosynostosis and identify a strong genetic influence on neurodevelopmental outcomes in affected children. These findings may eventually prove useful in determining which children with nonsyndromic craniosynostosis are most likely to benefit from surgical intervention. CLINICAL QUESTION/LEVEL OF EVIDENCE: Risk, III.

Unexpected role of SIX1 variants in craniosynostosis: expanding the phenotype of SIX1-related disorders.

BACKGROUND: Pathogenic heterozygous SIX1 variants (predominantly missense) occur in branchio-otic syndrome (BOS), but an association with craniosynostosis has not been reported. METHODS: We investigated probands with craniosynostosis of unknown cause using whole exome/genome (n=628) or RNA (n=386) sequencing, and performed targeted resequencing of SIX1 in 615 additional patients. Expression of SIX1 protein in embryonic cranial sutures was examined in the Six1nLacZ/+ reporter mouse. RESULTS: From 1629 unrelated cases with craniosynostosis we identified seven different SIX1 variants (three missense, including two de novo mutations, and four nonsense, one of which was also present in an affected twin). Compared with population data, enrichment of SIX1 loss-of-function variants was highly significant (p=0.00003). All individuals with craniosynostosis had sagittal suture fusion; additionally four had bilambdoid synostosis. Associated BOS features were often attenuated; some carrier relatives appeared non-penetrant. SIX1 is expressed in a layer basal to the calvaria, likely corresponding to the dura mater, and in the mid-sagittal mesenchyme. CONCLUSION: Craniosynostosis is associated with heterozygous SIX1 variants, with possible enrichment of loss-of-function variants compared with classical BOS. We recommend screening of SIX1 in craniosynostosis, particularly when sagittal±lambdoid synostosis and/or any BOS phenotypes are present. These findings highlight the role of SIX1 in cranial suture homeostasis.

Haploinsufficiency of SF3B2 causes craniofacial microsomia.

Craniofacial microsomia (CFM) is the second most common congenital facial anomaly, yet its genetic etiology remains unknown. We perform whole-exome or genome sequencing of 146 kindreds with sporadic (n = 138) or familial (n = 8) CFM, identifying a highly significant burden of loss of function variants in SF3B2 (P = 3.8 × 10-10), a component of the U2 small nuclear ribonucleoprotein complex, in probands. We describe twenty individuals from seven kindreds harboring de novo or transmitted haploinsufficient variants in SF3B2. Probands display mandibular hypoplasia, microtia, facial and preauricular tags, epibulbar dermoids, lateral oral clefts in addition to skeletal and cardiac abnormalities. Targeted morpholino knockdown of SF3B2 in Xenopus results in disruption of cranial neural crest precursor formation and subsequent craniofacial cartilage defects, supporting a link between spliceosome mutations and impaired neural crest development in congenital craniofacial disease. The results establish haploinsufficient variants in SF3B2 as the most prevalent genetic cause of CFM, explaining ~3% of sporadic and ~25% of familial cases.

Missense and truncating variants in CHD5 in a dominant neurodevelopmental disorder with intellectual disability, behavioral disturbances, and epilepsy.

Located in the critical 1p36 microdeletion region, the chromodomain helicase DNA-binding protein 5 (CHD5) gene encodes a subunit of the nucleosome remodeling and deacetylation (NuRD) complex required for neuronal development. Pathogenic variants in six of nine chromodomain (CHD) genes cause autosomal dominant neurodevelopmental disorders, while CHD5-related disorders are still unknown. Thanks to GeneMatcher and international collaborations, we assembled a cohort of 16 unrelated individuals harboring heterozygous CHD5 variants, all identified by exome sequencing. Twelve patients had de novo CHD5 variants, including ten missense and two splice site variants. Three familial cases had nonsense or missense variants segregating with speech delay, learning disabilities, and/or craniosynostosis. One patient carried a frameshift variant of unknown inheritance due to unavailability of the father. The most common clinical features included language deficits (81%), behavioral symptoms (69%), intellectual disability (64%), epilepsy (62%), and motor delay (56%). Epilepsy types were variable, with West syndrome observed in three patients, generalized tonic-clonic seizures in two, and other subtypes observed in one individual each. Our findings suggest that, in line with other CHD-related disorders, heterozygous CHD5 variants are associated with a variable neurodevelopmental syndrome that includes intellectual disability with speech delay, epilepsy, and behavioral problems as main features.

Exome sequencing implicates genetic disruption of prenatal neuro-gliogenesis in sporadic congenital hydrocephalus.

Congenital hydrocephalus (CH), characterized by enlarged brain ventricles, is considered a disease of excessive cerebrospinal fluid (CSF) accumulation and thereby treated with neurosurgical CSF diversion with high morbidity and failure rates. The poor neurodevelopmental outcomes and persistence of ventriculomegaly in some post-surgical patients highlight our limited knowledge of disease mechanisms. Through whole-exome sequencing of 381 patients (232 trios) with sporadic, neurosurgically treated CH, we found that damaging de novo mutations account for >17% of cases, with five different genes exhibiting a significant de novo mutation burden. In all, rare, damaging mutations with large effect contributed to ~22% of sporadic CH cases. Multiple CH genes are key regulators of neural stem cell biology and converge in human transcriptional networks and cell types pertinent for fetal neuro-gliogenesis. These data implicate genetic disruption of early brain development, not impaired CSF dynamics, as the primary pathomechanism of a significant number of patients with sporadic CH.

A genome-wide case-only test for the detection of digenic inheritance in human exomes.

Whole-exome sequencing (WES) has facilitated the discovery of genetic lesions underlying monogenic disorders. Incomplete penetrance and variable expressivity suggest a contribution of additional genetic lesions to clinical manifestations and outcome. Some monogenic disorders may therefore actually be digenic. However, only a few digenic disorders have been reported, all discovered by candidate gene approaches applied to at least one locus. We propose here a two-locus genome-wide test for detecting digenic inheritance in WES data. This approach uses the gene as the unit of analysis and tests all pairs of genes to detect pairwise gene × gene interactions underlying disease. It is a case-only method, which has several advantages over classic case-control tests, in particular by avoiding recruitment of controls. Our simulation studies based on real WES data identified two major sources of type I error inflation in this case-only test: linkage disequilibrium and population stratification. Both were corrected by specific procedures. Moreover, our case-only approach is more powerful than the corresponding case-control test for detecting digenic interactions in various population stratification scenarios. Finally, we confirmed the potential of our unbiased, genome-wide approach by successfully identifying a previously reported digenic lesion in patients with craniosynostosis. Our case-only test is a powerful and timely tool for detecting digenic inheritance in WES data from patients.